Infrared divergences of B meson exclusive decays to P-wave charmonia in QCD 

factorization and nonrelativistic QCD 
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In the framework of QCD factorization, we study the B meson exclusive decays B — > XcjK where 
the spin-triplet P-wave charmonium states Xcj{J = 0, 1, 2) are described by the color-singlet non- 
relativistic wave functions. We find that for these decays (except Xa) there are infrared divergences 
arising from nonfactorizable vertex corrections as well as logarithmic end-point singularities arising 
from nonfactorizable spectator interactions at leading-twist order. The infrared divergences due to 
vertex corrections will explicitly break down QCD factorization within the color-singlet model for 
charmonia. Unlike in the inclusive decays where the higher Fock states with color-octet cc pair and 
soft gluon can make contributions to remove the infrared divergences, their contributions can not be 
accommodated in the exclusive two body decays. As a result, the infrared divergences encountered 
in exclusive processes involving charmonia may raise a new question to the QCD factorization and 
NRQCD factorization in B exclusive decays. 
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B meson decays, especially exclusive nonleptonic de- 
cays, provide an opportunity to determine the parame- 
ters of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, 
to explore CP violation and to observe new physics ef- 
fects. However, although the underlying weak decay of 
the b quark is simple, quantitative understanding of non- 
leptonic B-meson decays is difficult due to the compli- 
cated strong-interaction effects. 

Exclusive £?-meson decays to final states containing 
charmonia are especially important due to three reasons. 
First, those decays e.g. B — > J/ij}K are regarded as 
the golden channels for CP studies due to their clean 
experimental signatures and straightforward theoretical 
interpretations. Secondly, those decays may provide use- 
ful information towards the understanding of color sup- 
pressed decays. Thirdly, those decays involve two energy 
scales, the beauty quark mass mb and charm quark mass 
m c , and therefore are more subtle in theoretical stud- 
ies. Since experimentally CLEO, BaBar and Belle Col- 
laborations have provided many measurements on the 
subject of B exclusive decays to charmonia [J such as 
J/ip,ip',rj c ,rj' c ,XcO, and Xci, theoretical studies on these 
issues become necessary. 

We will start with the QCD factorization approach 
proposed by Beneke et al.Q- It is argued that because 
the size of charmonium is small^l / a s m-ip) and its over- 
lap with the (B, K) system is negligible, the same QCD 
factorization method as for B —* 7T7t can be used for 
B — ► J/ipK decay. This small size argument for the 
applicability of QCD factorization to B exclusive de- 
cays to charmonia needs to be tested for specific decay 
channels. Indeed, explicit calculations for B — >■ J/ijjK 
within the QCD factorization approach 3] showed that 
the nonfactorizable vertex contribution is infrared safe 



and the spectator contribution is perturbatively calcu- 
lable at twist-2 order though the theoretical branching 
ratio is much smaller than the experimental data. Stud- 
ies on B — > rj c K decayQ also confirmed the qualitative 
applicability of the small size argument while the quan- 
titative underestimates compared with data were very 
similar to that encountered in B — > J/ijjK. However, 
further studies on B — > XcoK |5| challenged the applica- 
bility of QCD factorization for charmonia because of the 
appearance of non-vanishing infrared divergences in this 
decay. 

In previous studies for exclusive B decays to charmo- 
nia within the QCD factorization approach, the light- 
cone wave function description of the charmonia was 
adopted. This is fairly appropriate for the S-wave char- 
monia since the relative momentum q between charmed 
and anti-charmed quarks can be neglected in the lowest 
order approximation. However, for P-wave charmonia, 
q can not be neglected even in the leading-order. Al- 
though the light-cone wave function description of char- 
monia can provide us with some essential features of 
the problems involved in these decays, it can not give 
a complete analysis. The more appropriate method is to 
study charmonia within the non-relativistic bound-state 
picture or equivalently within the nonrelativistic QCD 
(NRQCD) framework 6J. In this letter, we report a com- 
plete analysis of the infrared unsafety encountered in B 
exclusive decays to P-wave charmonium states, includ- 
ing Xc2, which was not involved in previous studies, by 
using the color-singlet non-relativistic wave functions to 
describe the charmonia. 

To proceed, we first give a brief review of the QCD fac- 
torization approach. The general idea of the QCD factor- 
ization is that in the heavy quark limit mi> Aqcdj the 
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transition matrix elements of operators in the hadronic 
decay B — ► M X M 2 assumes a simple form at the leading 
order in l/mb|2]: 

(M 1 M 2 \O l \B)=F BAh {m 2 2 ) f dxT 1 (x)0m 2 (z) 

Jo 

fi 

d£ dx dy T n (£, x,y)<p B {&<j>M 1 {y)4>M 2 (x), (1) 

'o 

where Mi is the recoiled meson and M 2 is the emitted 
meson which is a light meson or a quarkonium. F BMl 
is the B — » Mi transition form factor, 4>m is the light- 
cone distribution amplitude and T 1,11 are perturbatively 
calculable hard scattering kernels. If we neglect strong 
interaction corrections, formula. reproduces the result 
of naive factorization. However, hard gluon exchange be- 
tween M 2 and BMi system implies a nontrivial convolu- 
tion of hard scattering kernels T 1 ' 11 with the distribution 
amplitude 4>m 2 - 

In the non-relativistic bound-state picture, instead of 
the light-cone amplitude distribution, charmonium can 
be described by the color-singlet non-relativistic wave 
function (the role of the higher Fock states with color- 
octet cc pair will be considered later on). Let p^ be 
the total 4-momentum of the charmonium and be the 
relative 4-momentum between c and c quarks. For P- 
wave charmonium, because the wave function at the ori- 
gin TZi(0)—0, which corresponds to the zeroth order in q 1 
we must expand the amplitude to first order in q. Thus 
we have (see, e.g., 0) 



where e a (L z ) is the polarization vector of an angular mo- 
mentum one system and the value of 7^(0) for charmonia 
can be found in e.g. Ref. 8] 

The spin projection operators Pss I (p, 0) and Pg S (p, 0) 
can be written as0 
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(7) 



[f(S z )ty + M) 7 a + 7 T {S z )(j/>+M)\ 



where we have made use of the non-relativistic approx- 
imation for the charmonium mass M ~ 2m. Here m is 
the charmed quark mass. 

In the calculation we need following polarization rela- 
tions for the 3 Pj states 

Y,t* a {L z y\S z ){lL z] \S z \W) = ±(- g ^ + P^-), 
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Y,^)e*\S z ){lL z ;lS z \U z ) - 
J2^ a (Lz)e* P (S z )(lL z ;lS z \2J z ) = e* af) (J z ), 



L Z S 



(8) 



LzSz 



where e\(J z ) is the usual spin-1 polarization vector and 
the polarization tensor e a ^(J z ) is that appropriate for 
a spin-2 system which is symmetric under the exchange 
a <->(3. 

The effective Hamiltonian relevant for B — > XcjK is 
written as0 



M(B ^Pj(cc)) (1L Z ; SS Z \JJ Z ) J ^Lq a H*^(y*VZC 1 Oi+ C 2 2 )- V»VZ £ C t O ( ) , (9) 

L *> s * v2 f=3 



x 5{q ~ l -^-WiM(<imO a (0)PssXp,0) + 0(O)P£ S >,O)] > ( 2 ) 

where 0{q) represents the rest of the decay matrix ele- 
ment and the spin projection operators Pss z (p, q) which 
is constructed in terms of quark and anti-quark spinors 
as 



p ss,(p,q)=\l— -<l,S2)u(~+q,si)(si;s 2 \SS z ), (3) 
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where Gf is the Fermi constant, Ci are the Wilson coef- 
ficients and V qiq2 are the CKM matrix elements. We do 
not include the effects of the electroweak penguin opera- 
tors since they are numerically small. Here the relevant 
operators Oi are given by 

01 = (s a bp)v~A ■ {c0C a )v~A, 

02 = (s a b a ) V - A ■ (cpCp)v-A, 

03(5) = (Sgb a )v-A ■ y^X<li3<lp)v-A(V+A)> (10) 



04(6) = (s a bp)v-A ■ ^(qpq^v- 



A(V+A) i 



o°(M do(q) \ 



dq a 

After q° is integrated out, the integral in Eq.© is pro- 
portional to the derivative of the P-wave wave function 
at the origin by 

/ ~0f qariMiq) = ^ a (L z )sj^ni{Q), (6) 



where a, (3 are color indices and the sum over q runs over 
u, d, s, c and b. Here {q\q 2 )v±A = 5i7m(1 =•= 7s)92- 

According to |2( all nonfactorizable corrections are due 
to Fig^ and other corrections are factorized into the 
physical form factors and meson wave functions. Tak- 
ing nonfactorizable corrections in Fig^into account, the 
decay amplitude for B —>XcjK(J = 0, 2) in QCD factor- 
ization is written compactly as 
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V cb V; s Ci-V tb V t * s (C 4 + C e ) xA, (11) 
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Xc2, 



32e;XKV^( 1 - z + ln2 ) 1 ,A 2 



m fe (l -z) 3 



ln( — *)+ finite terms, (16) 
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FIG. 1: Feynman diagrams for vertex and spectator correc- 
tions to B ^XcjK. 



where the coefficients A are given by 
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Here N c is the number of colors, C F = (N 2 - 1)/(2JV C ), 
and Fi is the B — > K form factor (see Ea. i|13|) below). 
The function // is calculated from the four vertex cor- 
rection diagrams (a, b, c, d) and /// is calculated from 
the two spectator correction diagrams (e, f) in FigQ] 
The form factors for B —>K are given as 



(K(p K )\s^b\B(p B )) 



(13) 



F^Kpb+Pk), + [F (p 2 ) - Fxip 2 )] ^ m K 

where p = p B —px is the momentum of charmonium with 
p 2 = M 2 , and m B , are respectively the masses of B, 
K mesons. We will neglect the kaon mass for simplicity. 
We can use the ratio between these two form factors as 
F (p 2 )/Fi(p 2 ) =l—p 2 / m 2 B . So we need only one of the 
two form factors, say, F\ to describe the decay amplitude 
just like that in Ea. l(T2|) . 

Since our main task in this letter is to explore the in- 
frared unsafety due to the vertex corrections, we only 
give the explicit results below for the infrared divergence 
terms in // which are regularized within the gluon mass 
scheme. We also present the leading- twist results for ///. 
Here // and fjj correspond respectively to the T 1 ,T IT 
terms in Eq.Q. 

For 3 Pi charmonium state \ci, there is no infrared di- 
vergence in B — > Xc\K decay, and we will give the rea- 
soning shortly (results of this decay within the light-cone 
wave function description can be found in Ref. :5| ) . 

For 3 Po charmonium state \c0, 



8rribz(l — z + lnz) A 2 

// = ~ — ln (-2 ) 



(1 



finite terms, (14) 



Here z — M 2 /m 2 B « 4m 2 jm\, £ is the momentum frac- 
tion of the spectator quark in the B meson and y = (1— y) 
is the momentum fraction of the s quark inside the K 
meson, 4>b and 4>k are the light-cone wave functions for 
the B and K mesons respectively 2] . Here A is the gluon 
mass introduced to regularize the infrared divergences in 
the vertex corrections. We have simplified the results for 
/// by noting that £ ~ 0(Aq C dA™&) ~> °- 

The integrals for /// in Eas. lfT3|l . ((T7|) will give loga- 
rithmic divergences when the asymptotic form 4>k(_u) = 
6y(l — y) for the kaon twist-2 light-cone distribution am- 
plitude as well as the same parameterization for the £ 
integration as that in Ref.|2| are used. It is also worth- 
while to emphasize that the infrared singularities in fx 
are more serious than the end-point singularities in ///. 
This is because, when one considers the effects of parton 
transverse degrees of freedom for B and K mesons, end- 
point singularities arising from spectator interactions can 
be regularized |10| and logarithmic divergences can then 
be removed. However, the infrared divergences due to 
nonfactorizable vertex corrections still exist. 

Since two large scales are involved, we can have sev- 
eral different choices of the heavy quark limit. In the 
above calculations, we have chosen the heavy quark limit 
as Wfc — >oo with m c /mb fixed. In this limit, the infrared 
divergences in // as well as the logarithmic divergences in 
fn, will break down QCD factorization. Another choice 
of the heavy quark limit is that mj, — > oo with m c fixed. 
Then all the divergences mentioned above are power cor- 
rections and should be dropped out, so QCD factoriza- 
tion still holds in this limit. Physically, the latter case is 
equivalent to the limit of zero charm quark mass in which 
charmonium is regarded as a light meson. Obviously, the 
first choice of the heavy quark limit with m c /mb fixed 
is more interesting in phenomenological analysis, and is 
also usually used in theoretical studies 0, EL where it 
is expected that QCD factorization should apply to B 
meson exclusive decays into charmonium in the limit 
m c —>■ oo with corrections of order Aqcd /(m c a s ) ~ 1. 
Our result shows that this expectation holds for decays 
to S-wave charmonia but not for P-wave charmonia we 
studied above where the vertex infrared divergence will 
break down factorization at order of AQCD/(w c a s ). 
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It is worthwhile to discuss the origin of the infrared 
divergences and the differences between S-wave and P- 
wave charmonia as well as the differences among P-wave 
charmonium states. We write momenta for cc as p c — 
p/2 + q, pc = p/2 — q. The gluon coupling to the cc pair 
showing in Figgis given by 



Ju = 



"/u($c + ft + m c )T 
(p c + k) 2 - m 2 



^5 + $ - m c )^ v 
(pc + k) 2 — m 2 



(18) 



where T denotes for the weak decay vertex. When k 
is soft, the coupling in the nonrelativistic expansion in 
terms of q ■ k/p ■ k can be simplified to 



</„ ~ 4r 



q v (q ■ k)p v 



p ■ k (p ■ k) 2 



(19) 



where we used the on-shell conditions for c and c quarks. 
For S-wave charmonium states such as J/tp, since the 
wave functions depend ong 2 , q in Eq. (|19f) makes no con- 
tribution to the lowest order. So there is no infrared 
divergence for S-wave charmonium states. For Xd, be- 
cause the orbital and spin momenta couplings are in an 
anti-symmetric form (see Eq.(jHJ), when combined with 
the B — > K form factors the infrared divergences aris- 
ing from Eq. (|19fl are totally cancelled out (note that this 
cancellation only accidentally holds because of the spe- 
cific form of B — > K form factors) . However, for XcO and 
X C 2, because the orbital and spin momenta couplings are 
in a symmetric form, infrared divergences arising from 
Eci. (|19fl are not cancelled out. This qualitative argument 
is confirmed by our explicit calculations performed above. 

It is well known that there are infrared divergences in 
the inclusive decay and production of P-wave charmonia, 
which are related to the problems encountered here. In 
the inclusive processes, the infrared divergences in the 
color-singlet P-wave cc state can be removed by including 
contributions from the higher Fock states with color-octet 
cc pair (say in S-wave) and soft gluon within the NRQCD 



factorization framework (e.g. for B decay see [T^l, for an- 
nihilation hadronic decay see [l3L Il"H , and for decay and 
production see 0, Hfil ll^|). However, to the best of 
our understanding, the color-octet cc pair with dynami- 
cal soft gluon can make contributions to the multi-body 
but not two-body exclusive decays. (Note that the soft 
gluon emitted by the color-octet cc and subsequently re- 
absorbed by the spectator light quark has already been 
considered in the nonfactorizable spectator contribution 
discussed above (see diagrams (e,f) in Fignj, where the 
soft gluon is related to the end-point singularities of B 
and K light-cone wave functions.) As a result, it seems to 
us that the infrared divergences encountered in exclusive 
processes involving charmonia may raise a new question 
to the QCD factorization and NRQCD factorization in B 
exclusive decays. Further studies are needed to seek the 
solution to remove the infrared divergences. 

In summary, we have studied the exclusive two-body 
decays of B meson into the spin-triplet P-wave charmo- 
nium states Xcj{J =0, 1, 2) and kaon within QCD factor- 
ization and NRQCD by adopting the color-singlet non- 
relativistic wave function for charmonium. We find that 
for these decays (except B — > \ c iK), there are infrared 
divergences arising from nonfactorizable vertex correc- 
tions, which can not be removed by including contribu- 
tions from the color-octet cc with soft gluon Fock states, 
and therefore will break down QCD factorization. New 
considerations should be introduced to describe B meson 
exclusive decays to charmonium states. 
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